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ABSTRACT 

ON 

A formula for off-axis emissions from the beamed afterglow of GRBs in a synchrotron 
model is derived. The formula is applied to an inhomogeneous circumstellar matter 
obeying n oc r~ 2 suggested for SN1998bw/GRB980425. The slow decline rate as well 
\q . as the X-ray flux similar to the observed values are obtained for appropriate choice 

of the parameters so that the association of GRB980425 with SN1998bw is further 
strengthened. 

(ApJ Letter in press) 
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have been reported ( Pian et al. 19991 ), The NFI observation identified an X-Ray source 



Subject headings: gamma rays: bursts — stars: supernovae 



1. Introduction 



Recently the X-ray light curve of GRB980425 and follow up data obtained by the BeppoSAX 



(SI) position-ally coincident with the unusual Type Ic supernova SN1998bw ( Galama et al 



1998| , [Kulkarni et al. 1998 ). Note here that the position of SI is revised and the other source S2 



is ~ 4 away from SN1998bw. The X-ray flux of SI was 3 x 10 13 erg s 1 cm 2 in April 26 and 
declined a factor of two in six months. If this is an X-ray afterglow of GRB980425, the decline 
rate (~ t~ ' 2 ) is extremely slower than usual, which needs an interpretation. If GRB980425 is 
associated with SN1998bw, the distance is ~ 40Mpc ( Palama et al. 1998 ) and the isotropic 



gamma ray total energy is unusually low ~ 10 48 erg compared with > 10 51 erg in other GRBs for 
which a redshift measurement is available. 

Such a low gamma ray luminosity might be interpreted if the beam direction of GRB980425 
is different from the line of sight ( [Nakamura 1998| , Eichler fc Levinson 1998 ). In such an geometry 



we see only scattered gamma rays so that the luminosity is low and the energy of gamma rays is 
not high < 500keV (Nakamura 1995, Eichler fc Levinson 1998] ). This also interprets the fact that 



GRB980425 belongs to so called Non High Energy class of GRBs flPendelton et al. 1997| ). It is 



shown that if the angle between the beam and the line of sight is ~ 30°, the luminosity can be 
~ 10 -3 smaller than usual (Eichler & Levinson 1998). Such a beam model is relevant now since an 



evidence for the beaming has recently been suggested from the rapid decline of optical afterglow 
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of GRB990123 after day 2 flKulkarni et al. 1999|) and the rapid decline of GRB980519 Q Halpern 
Kemp, Piran fc Bcrshady 1999| ), 



An angle of ~ 30° has also been suggested to interpret the linear optical polarization (~0.5% 
in June 20 ~ 23; Kay, Halpern fc Leighly, 1998] ) of SN1998bw . If the supernova ejecta is an 
oblate or a prolate shape with an aspect ratio 1 : 2 ~ 3 and the line of sight is ~ 30° away from 
the symmetry axis, the observed linear polarization degree might be explained ( Hoeflich, Wheeler 



fc Wang 1998|) . 



From the radio observation of SN1998bw it is suggested that the density distribution of the 
circumstellar matter follows n = nQr 2 /r 2 with nor 2 , ~ 10 34 /cm, which is consistent with the 
constant mass loss rate in the progenitor ( Li fc Chevalier 1999[ ). It is also claimed that the energy 
of the order of 10 50 erg should be injected twice with the injection velocity of the order of 0.5c, 
which strengthens the link between SN1998bw and GRB980425. 

In this Letter we try to interpret the slowly declining X-ray afterglow of GRB980425 in the 
beam model of GRBs assuming the line of sight is ~ 30° away from the axis of the beam. We 



adopt the density distribution of the circumstellar matter suggested by Li fc Chevalier 199£ . 



We will obtain the decline rate and the X-ray flux similar to the observed values for appropriate 
choice of the parameters so that the association of GRB980425 with SN1998bw will be further 
strengthened. 



2. Off-axis emission from the beamed afterglow 

I ~l . _. 

Granot, Piran & Sari 1998 as well as Woods &: Loeb 1999| derived a general formula 



to compute the off-axis emission from beamed GRBs. Here we adopt their formulations and 
notations. Let us use a spherical coordinate system r = (r,9,<j)) where the coordinate are measured 
in the lab frame; let the 9 = axis (z-axis) points to the detector and r = be the central engine. 
Let also D be the distance to the source and a = rsin0/D be the angle that a given ray makes 
with the normal to the detector. Then the observed flux is given by 

F V {T) = — / d(p / a da — , , (1 

7/3 Jo Jo Jwy(i-p) v u {1 - n 2 } 6 / 2 

where [i = (1 — u'/ r yu)/f3. a m , T and j' u , are the maximum value of a, the arrival time of a photon 
at the detector and the rest frame emissivity measured in erg s _1 cm -3 Hz -1 sr" 1 , respectively. 
Note here that ' means the physical quantity in the rest frame. 



We adopt the standard fire ball model( Piran 199"8| ). As an emissivity we extend a simple 



synchrotron model( |Sari, Piran fc Narayan 1998| ) used in a homogeneous ambient gas to the 
inhomogeneous circumstellar matter following n = no(r / ro)~ d . d=2 corresponds to the constant 
mass loss from the progenitor of GRBs. This kind of a model has ever been studied in a different 
context ( pVleszaros, Rees fc Wijers 1998j ). 
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We assume that electrons are accelerated in the shock to a power-law distribution of 
Lorenz factor 7 e with a minimum Lorenz factor N(j e ) oc 7~ p d7 e ,7 e > 7 m . 7 m = G7 
where 7 is the Lorenz factor of the shocked fluid and G = e e (p — 2)/(p — l)m p /m e with e e 
being the efficiency of the acceleration. Since p = 2 ~ 2.5 is suggested from observations, 
G = 60(e e /0.1)(3(p — 2)/(p — 1)). The magnetic field strength is given as B = Boj(r /ro)~ d l 2 
where Bq = v /327rm p noe_B. £b measures the ratio of the magnetic field energy to the total 
thermal energy. The rest frame radiation power (P 1 ) and the characteristic synchrotron frequency 
(z/(7 e )) from a randomly oriented electron with Lorenz factor 7 e > 1 in a magnetic field B 
are given by P' = PoJ 2 Je( r / r o)~ d an d ^'(7e) = v oll1 (r/ro)~ d / 2 where Pq = o~tB 2 c/6it and 
z/ = eB /2irm e c = 3.5 x 10 5 Hz(n /lcm- 3 ) 1 /2( eB /o.l)V2. 

The peak spectral power flux does not depend on j e and occurs at f'{j e ) and is 
given by P '' u , max = P'/u' = P®y /vo(r /r$) ~ d l 2 . The critical gamma factor 7 C is defined by 
j c m e c 2 = P ! {^ c )t = P' (j c )r / (-yc) and is given by 7 C = F '/"f(r /ro) d ~ l where F = m e (? / 'PqTq. Then 
the energy spectra depend on two frequencies defined by v' m = 1^(7™) = uqG 2 ^ 3 (r /ro)~ d / 2 and 
' c = v '{ lc ) = vqF 2 ^ 1 ^/^)^/ 2 - 2 ( Pari, Piran fc Narayan 19981 ) . 



There are two different cases in the energy spectra depending on whether v' m > v' c (fast 
cooling) or v' m < u' c (slow cooling); ( pari, Piran &: Narayan 1998j ). In both cases the emissivity j' v , 
is expressed in the form as 

ji' = K', m ^(-r d f^), (2) 

47T Tq 

where f(v') is given for fast cooling case as 

( (v'/v'J 1 / 3 v' c > v' 

fW) = WK)~ 1/2 V'm>v'> V'c 

{ KJ<)- 1/2 (v'/v' m )- p/2 V > u' m , 
while for slow cooling case we have a smilar formula( |5ari, Piran &: Narayan 1998 ). 



Let A9 and 9 V be the beaming half-angle and the angle between the direction to the detector 
and the axis of the emission cone. In case of 9 V = 0, Eq. (1) is evaluated as 

K{T) = 1 = r~" s&wm- <3) 

For 7 ^> 1 and A9 < 1 , the integral is expressed as 

n 2 Ae 2 js 

^lP^ ma J(vis/2). (4) 

If 7A0 ^> 1 , I does not depend on A9. This is expected since in this case the detector detects the 
radiation only from the half-angle of 7" 1 and can not observe the edge of the beam. For d = 0, 
Eq. (4) gives the same results as in Sari, Piran fc Narayan 1998| . \_ 



In reality there is 10% or so difference depending on the parameter p. However such a difference is consistent in 



a simple model of emissivity in Sari, Piran & Narayan 1998 
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Now let us consider another extreme case of 1 > 9 V > AO > 7 1 . In this case the line of sight 
from r = does not pass through the beam. Then Eq. (1) is evaluated as 

FJT) _ ±^ { L f ^n h j _ (5) 

\1txD a vq r J V 7 (i_ p cos(6»„-Ae)) 1/ 

where cos ^„ = (cos A9 — cos 6* cos flu) /sin 9 V sin 

Eq.(5) is a general formula for off-axis emission from beamed GRBs in a synchrotron model. 
In the next section we shall apply the formula to GRB980425. 



3. Application to GRB980425 



Now let us consider d=2 case which is suggested from SN1998bw ( Li fe Chevalier 1999 ). 
Then Eq. (5) is rewritten as 

F ^ = 1^7 /; 1 = / 2^f(vjs). (6) 

11ttU z Uq J(1-Pcos{O v -A0)) s 

The position of the shock is related to the arrival time T as r = cT/(l — /3cosfl„) ~ 
1.9 x 10 16 cm(fly/30°) _2 Trf where is the arrival time in days. 7 is given for the adiabatic case by 

, ^ a/2 = go/ #6 -/^cosfl^ 1/2 n 1/2 r x 1/2 

7 l m p c 2 n r^27rA^ ; l 10 51 er g ; 1 A# 2 j 4cm-3 J l 10 17 cm J d ' ^ ' 

where is the total energy in the beam. 

The two characteristic energy hv m and hv c at the detector are given by 

hu m = 2x 10- 2 eV(g) 2 ( lTl^ )-i ( ^ )2( ^ ) -i ( ^ ) i/2 ( ^ ) i/2 (8) 
60 0.1 20 ro 0.1 1cm 13 

fo/ c = 8.7 x io-3 e V( 1 -^ g " )-H^)- 2 (-)(^)- 3 / 2 ( T ^)- 3 / 2 (^)- 2 0) 

0.1 20 ro 0.1 1cm 6 10 i( cm 

Now consider the observed radiation in optical and X-ray bands. I is evaluated for 9 V ~ 30° 
and fast cooling case as 

/ = 2 ^+p/ 2 4^(—) 2 G p 'H L ) l ~ p/2 Y' 2 n-) p/2 - (10) 

9 V +P flu r v 

Then the energy flux in a band of v& < v < v u is given by 

F b = r F v du = 1.74 x 10- 12 erg s" 1 cm^O,-?-^)*'/ 2 - 1 - (^f/ 2 " 1 ) 
Ju d P - 2 v d v u 

( r \4-p ( n \2-p/2 ( D \-2 ( °v N-4/ A6l N4-p r E b sp/l-lrpl-p ( , \ 

l 10 17 cm J 4cm-3 ; M0Mpc j { 30 o) { 9 V ' l 10 51 erg j d ' { } 
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where C p = 2 3 (P" 2 )/ 2 12140( p " 2 ) (G'/60) p " 1 Let us compute F b in X-ray band (y d = 1.6keV, 
v u = lOkeV) and V band (u d = 5x 10 14 Hz, v u = 6 x 10 14 Hz) for p = 2.2 and p = 2.5. 

For p = 2.2 

F b = 6.5 x 10~ 13 erg s" 1 cm~ 2 (1.4 x 10~ 13 erg s" 1 cm" 2 ; V band)( 



10 17 crn 

lcm 3 40Mpc 30 0.1 0.1 t> v 10 Oi erg 

In this case the declining rate is almost the same as the observed value while the absolute value 
of the X-ray flux is similar to the observed one and can be adjusted by choosing the appropriate 
values of various parameters such as 6 V , A8, e e , e#, r , n and E b . 

For p=2.5 

\1.5 



F b = 8.0 x 10" 13 erg s" 1 cm" 2 (3.0 x 10" 13 erg s" 1 cm" 2 ; V band)(— ^- , 

10 i 'cm 

/_^0 \0.875/_^_\-2 ( ' 9v n-4/ e e s 1.5 / e B \0.125 ( ^0 \1.5 /_^£> ^0.25^-0.5 /io\ 

4cm- 3j l 40Mpc J l 30° j l 0.1 j l 0J j 1 B v ] [ W^erg> ±d ' [L6) 

In this case, the amplitude of the X-ray flux can be similar to the observed value but the decline 
rate is somewhat too rapid. 

For comparison, we shall compute the energy flux for an aligned case (6 V = 0). Since in this 
case the position of the shock is related to the arrival time of the radiation by r = 2j 2 cT, 7 is 
given by 

ry = 1 3 ( Eb )l/4 ( ^)-l/2 ( n n-1/4, ^0 )-l/2 T -l/4 (14) 

7 l 10 51 erg j l 15° ; { lcm-^ l 10 17 cm J d 1 ' 



F b is given for p = 2.2 as 

F b = 7.8 x 10~ n erg s" 1 cm" 2 (1.7 x 10" n erg s" 1 cm" 2 ; V band)( : 



r o ^0.5 



10 17 cm' 

1 n ^0.3/ D \-2 / AO \-1.6 ( £ e n 1.2 ^ € B xQ.l/ Eb ^0,8^-0.9 r\ 

4cm- 3 J l 40Mpc ; K 15°' { 0.l' l 10 51 erg J d ' K } 



while for p = 2.5 



r x0.5 



F b = 2.7 x 10" n erg s" 1 cm" 2 (1.7 x 10" n erg s" 1 cm" 2 ; V band)(- - 

10 i 'cm 

i n ^0.375/ D \-2( ^ n — 1.75 / e e \1.5/ e B \Q.125/ ^ b nQ. 875^-1. 125 

4cm- 3 J l 40Mpc J l 15° J { 0.l ) { 0.1 } 4o 51 erg j d ' 1 ' 

Let us apply Eq. (16) to GRB990123. There was a break in the optical light curve at T d = 2. 
This may be due to the widening of the beaming half-angle ( Rhoads 1997 , Kulkarni et al. 1999| ). 



From the condition of 7 = A9 at T d = 2, A9 is expressed by the other parameters as 
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Adopting D =10Gpc as the luminosity distance to GRB990123 (z = 1.6 ; Bloom et al. 1999a| ) we 
have 

F b = 1.8 x 10" 14 erg s" 1 em _2 (l.l x 10" 14 erg s" 1 cm" 2 ; V band)(^J — )" L125 

10 i 'cm 

( n \-0-5/ D \-2f e e \1.5f e B \0.125/_^> \1. 75^-1. 125 / 1S \ 

( i^ ) ( ioG P ^ ) ( oT } ( oT ) ( io^i ) Td ■ (18) 

For example, if Ei ~ 10 52 erg with the other parameters unchanged, the decline rate and amplitude 
of the flux in X-ray and the optical bands are similar to the observed values. However even if 
Ef, ~ 10 51 erg, by appropriate choice of the other parameters such as r , ra , e e and e#, we may also 
have the observed amplitude. 



4. Discussion 

Since a possible association of GRBs with Type Ib/Ic supernovae is suggested only for 
GRB980425, we had better wait for more events to confirm the conjecture. However recently 
another possible association of GRB with the supernova is suggested for GRB980326. The decline 
of the afterglow was leveled off at ~20 days after the GRB, which was considered as the host 
galaxy. Keck observed the host galaxy nine months after the burst again but could not find the 
host galaxy. A possible explanation for this peculiar event is that the apparent level off at ~20 
days is due to the peak of supernova like SN1998bw([Bl oom et al. 1999b| ). If this is the case, at 
least two examples of association of GRBs with supernovae exist. 

For GRB980425 in our model the break in the light curve when 7 = 8~ l should occur at 

Td = 5 44(^^)(^) 2 (^ T )- 1 (^^)- 2 . (19) 
V 10 51 erg n 30° y V lcm~ 3; V 10 17 crn V ; 

The break time does not depend on A6 but on Ef,,6 v ,no and ro- Before the break time, Eq.(12) 
and (13) can be used so that the optical luminosity is much lower than the luminosity of SN1998bw 
for Trf < 180 ( Iwamoto et al. 1998| ). However the decline rate is so low that the luminosity from 



the afterglow might overcome the luminosity of the supernova and the host galaxy before Td ~ 544. 
Therefore it is important to observe SN1998bw up to ~ 2 years after the supernova event. 

In conclusion, if the line of sight is ~ 30° away from the axis of the beam for GRB980425, the 
decline rate and the X-ray flux similar to the observed values are obtained for appropriate choice 
of the parameters. This strengthens the link between GRB980425 and SN1998bw. 
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